Molecular Systems Biology 9; Article number 632; doi:10.1038/msb.2012.65 ^ moleCUlar 

Citation: Molecular Systems Biology 9:632 ClJ W SyStemS . 

© 2013 EMBO and Macmillan Publishers Limited All rights reserved 1744-4292/13 open transparent blolOgy 

www.molecularsystemsbiology.com access process 

An Oct4-Sall4-Nanog network controls developmental 
progression in the pre-implantation mouse embryo 

Meng How Tan^ '^^ Kin Fai Au^'^ Denise E Leong\ Kira FoygeP '^ Wing H Wong^'^ * and l\/lylene Wl\/I Yao^'^ * 

^ Department of Obstetrics and Gynecology, Stanford University School of Medicine, Stanford, CA, USA, ^ Department of Statistics, School of Humanities and 
Sciences, Stanford University, Stanford, CA, USA and ^ Department of Health Research and Policy, Stanford University School of Medicine, Stanford, CA, USA 
"^Present address: Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA 

^Present address: Department of Radiology, Molecular Imaging Program, Stanford University School of Medicine, Stanford, CA 94305, USA 
^Present address: Univfy Inc., Los Altos, CA 94022, USA 
^These authors contributed equally to this work 

* Corresponding authors. WH Wong, Department of Statistics, School of Humanities and Sciences, Stanford University, 390 Serra Mall, Stanford, CA 94305, USA. 
Tel.: +1 650 725 2915; Fax: +1 650 725 8977; E-mail: whwong@stanford.edu or MWM Yao, Department of Obstetrics and Gynecology, Stanford University School of 
Medicine, Stanford, CA 94305, USA. Tel.: +1 650 799 8003; Fax: +1 650 961 1377; E-mail: mylene.yao@gmail.com 

Received 28.7.11; accepted 30.11.12 



Landmark events occur in a coordinated manner during pre-implantation development of the 
mammalian embryo, yet the regulatory network that orchestrates these events remains largely 
unknown. Here, we present the first systematic investigation of the network in pre-implantation 
mouse embryos using morpholino-mediated gene knockdowns of key embryonic stem cell (ESC) 
factors followed by detailed transcriptome analysis of pooled embryos, single embryos, and 
individual blastomeres. We delineated the regulons of Oct4, Sall4, and Nanog and identified a set of 
metabolism- and transport-related genes that were controlled by these transcription factors in 
embryos but not in ESCs. Strikingly, the knockdown embryos arrested at a range of developmental 
stages. We provided evidence that the DNA methyltransf erase DnmtSb has a role in determining the 
extent to which a knockdown embryo can develop. We further showed that the feed-forward loop 
comprising DnmtSb, the pluripotency factors, and the miR-290-295 cluster exemplifies a network 
motif that buffers embryos against gene expression noise. Our findings indicate that Oct4, Sall4, and 
Nanog form a robust and integrated network to govern mammalian pre-implantation development. 
Molecular Systems Biology 9: 632; published online 8 January 2013; doi:10.1038/msb.2012.65 
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Introduction 

Many critical cellular events occur within a few days during the 
pre-implantation development of a mammalian embryo. First, 
the paternal and maternal genomes are reprogrammed to a 
totipotent state. Second, the embryonic genome is activated 
after a period of reliance on maternal factors. Third, polarity is 
established when the inner and outer cells of a multicell 
embryo are biased toward different cell fates. Finally, two 
important cell fate decisions occur before implantation. In the 
first cell fate decision, the outside cells develop into extra- 
embryonic trophectoderm, while cells that are inside the 
embryo retain pluripotency and form the inner cell mass (ICM) . 
In the second cell fate decision, some of the ICM cells further 
differentiate into the primitive endoderm (Zernicka-Goetz etaU 
2009). Given their complexities, all of these early cellular 
events have to be carefully orchestrated by a gene regulatory 
network to ensure the proper development of the embryo. 

The gene regulatory network of embryonic stem cells (ESCs) 
has been extensively investigated in recent years. ESCs, which 
are derived from the ICM of mammahan embryos, are 
pluripotent like cells in the ICM but they can also self-renew 
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indefinitely in culture. The transcription factor Oct4, also 
known as PouSfl, has been shown to have a key role in the 
maintenance of an undifferentiated state in both human and 
mouse ESCs. Chromatin immunoprecipitation (chIP) assays in 
ESCs reveal that Oct4 binds to more than a thousand loci in the 
genome, suggesting that it may regulate many genes as a 
guardian of pluripotency (Boyer et aU 2005; Loh et aU 2006; 
Chen et aU 2008a} . In addition, other key transcription factors 
that control the self-renewing and pluripotent properties of 
ESCs include Sall4 and Nanog. Previous studies indicate that 
Oct4, Sall4, and Nanog share a close functional relationship. 
Co-immunoprecipitation experiments show that they interact 
with one another in vivo (Wu et aU 2006; Liang et aU 2008; 
Yang et aU 2008). Furthermore, the three transcription factors 
co-occupy the promoters of many target genes (Boyer et aZ, 
2005; Wu et aU 2006; Yang et aU 2008), indicating that they 
form an integrated network in ESCs. 

Compared with ESCs, much less is known about the gene 
regulatory network in early mammalian embryos largely due 
to technical challenges. In particular, limited number of 
embryos or cells has so far precluded extensive studies of 
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any regulatory network operating during pre-implantation 
development. For example, chIP experiments to map tran- 
scription factor binding sites in pre-implantation embryos are 
currently unfeasible. Fortunately, advances in gene expression 
profiling are now allowing studies into small number of cells or 
even single cells. Several laboratories used microarrays to 
profile the transcriptome during pre-implantation develop- 
ment of mouse embryos (Hamatani et aZ, 2004; Wang et aU 
2004; Zeng et al 2004; Xie et al 2010). Expression analysis of 
48 genes in individual cells from the one-cell zygote to the 
blastocyst identified markers of outer and inner cells (Guo 
et al 2010} . In addition, next-generation sequencing technol- 
ogy was used to investigate the pre-implantation transcrip- 
tome in single mouse blastomeres (Tang et al 2011). While 
these previous studies were important for understanding early 
mammalian development, they relied primarily on transcrip- 
tome profiling of wild-type embryos and did not perform 
perturbation experiments to dissect interactions between the 
different genes. 

Here, we show that the key pluripotency factors in ESCs also 
form critical nodes in the regulatory network governing pre- 
implantation development of the mammahan embryo. We 
used gene-specific knockdowns followed by genome-wide 
microarray analysis as well as extensive single-embryo and 
single-cell expression profiling to deconstruct the Oct4-Sall4- 
Nanog transcriptional network in pre-implantation embryos. 
Our study yielded novel and unexpected insights into the 
regulatory mechanisms that control the earhest stages of 
mammahan development. 

Results 

Oct4, Sall4, and Nanog are essential for pre- 
implantation development of the mouse embryo 

Since a network of transcription factors has been shown to 
regulate pluripotency in ESCs (Zhou et al 2007; Chen et al 
2008a}, we hypothesize that some of these factors may also 
have an integral role in controlling developmental progression 
of the pre-implantation embryo. We previously found that 
Oct4 may play an important role at the maternal-zygotic 
transition (Foygel et al 2008} . Here, we sought to undertake a 
more comprehensive study of the functions of key ESC factors 
in pre-implantation embryos. We first measured the abun- 
dance of several pluripotency factors in individual one-cell 
zygotes using a microfluidic real-time PCR system known as 
Biomark (Fluidigm}. While the Kriippel-like factor KlfS was 
absent, we detected the transcripts of Oct4, Sall4, and Nanog in 



multiple one-cell embryos (Supplementary Figure SI}. Hence, 
we chose to focus on the three transcription factors in this study. 

To determine if Oct4, Sall4, or Nanog played an important 
role during pre-implantation development of the mammalian 
embryo, we depleted each of the proteins by injecting 
translation-blocking morpholinos (MOs} (Supplementary 
Figure S2}, which targeted both maternal and embryonic 
transcripts, into one-cell fertilized zygotes and then cultured 
the injected embryos together with control embryos in serum 
supplemented Human Tubal Fluid (SAGE In- Vitro Fertilization 
Inc.}. For each MO injection, we confirmed that the protein 
was absent by immuno cytochemistry (Figure 1 A} . After 4 days 
of in vitro culture, 97.7% of uninjected embryos and 85.1 % of 
control MO-injected embryos developed at least until the 
morula stage, but surprisingly, <25% of embryos that were 
depleted of any of the three pluripotency factors compacted. In 
addition, 40-60% of embryos injected with a MO targeting 
Oct4, Sall4, or Nanog arrested by the 4-cell stage after 4 days of 
in vitro culture, compared with <10% for the uninjected or 
control MO-injected embryos (P<0.05, Student's Mest}. To 
verify that our findings were not due to off-target effects, we 
performed rescue experiments for each pluripotency factor 
(Figure IB and C}. The percentage of embryos that arrested by 
the 4-cell stage decreased by 25.3 % when we co-injected Oct4 
MO with Oct4 mRNA, while the percentage decrease for Sall4 
or Nanog was 18.3 or 24.7% respectively (P<0.05, Student's 
t-iesi) . These and other additional rescue data (Supplementary 
Results and Supplementary Figures S3 and S4} indicated that 
the MOs were specific for their intended targets and that our 
observed phenotypes were a direct consequence of knocking 
down the pluripotency factors. Taken together, our experi- 
ments indicate that Oct4, Sall4, and Nanog are essential for 
pre-implantation development. Interestingly, these results 
contradicted pubhshed knockout mouse models, where 
embryos that were homozygous null for Oct4, Sall4, or Nanog 
appeared to develop in vivo at least until the blastocyst stage, 
although they were still embryonic lethal (Nichols et al 1998; 
Mitsui et al 2003; Sakaki-Yumoto et al 2006}. A likely reason 
for the apparent discrepancy is that maternal transcripts are 
not eliminated in knockout mouse models (see Supplementary 
information for more discussion} . 



Global transcriptome analysis identifies genes 
dependent on Oct4, Sall4, and Nanog 

To understand the molecular mechanisms that might cause the 
embryos to arrest during knockdown of a pluripotency factor. 



Figure 1 MO-mediated knockdown of Oct4, Sall4, or Nanog in pre-implantation mouse embryos was effective and specific and resulted in developmental arrest. 

(A) Confocal microscopy images of embryos stained with a-0ct4 (left panel), a-Sall4 (middle panel), or a-Nanog (right panel) antibodies. The embryos were stained for 
the protein-of-interest 2 days after microinjection. Oct4 and Nanog proteins were completely undetectable by immunostaining in 100% of the embryos injected with the 
corresponding MO, while the uninjected embryos and the control MO-injected embryos were strongly stained. In addition, while a small amount of Sall4 protein was 
detected in embryos injected with a Sall4 MO, the protein level in these embryos was significantly less than that in all the control embryos. (Scale bar = 45|im). 

(B) Representative images of pre-implantation embryos under various experimental conditions. For every set of rescue experiments, there were four conditions: (i) no 
injection, (ii) injection with a control MO, (ill) injection with the experimental MO alone, and (iv) co-injection with the experimental MO and the corresponding mRNA that 
was in vitro transcribed. (C) Percentages of embryos that arrested by the 4-cell stage during 4 days of in vitro culture. Left panel (green bars): compared with embryos 
injected with the Oct4 MO only, we observed that co-injecting 0.6 mM Oct4 MO with 0.20-39 ng/|il Oct4 mRNA resulted in a smaller percentage of embryos arresting by 
the 4-cell stage (9 independent experiments) (P<0.05, Student's f-test). Middle panel (blue bars): compared with embryos injected with the Sall4 MO alone, we found 
that co-injecting 0.6 mM Sall4 MO with 7-146 ng/|il Sall4 mRNA lowered the rate of developmental arrest by the 4-cell stage (8 independent experiments) (P<0.05, 
Student's f-test). Right panel (red bars): compared with embryos injected with the Nanog MO alone, we discovered that co-injecting 0.6 mM Nanog MO with 5 ng/|il 
Nanog mRNA also resulted in partial rescue of the phenotype (7 independent experiments) (P<0.05, Student's Mest). 
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we investigated transcriptome changes in pooled embryos 
using Affymetrix microarrays at ~20 and 44 h after MO 
injection (Supplementary Files S1-S5). At 20 h, uninjected 
control embryos are at the 2-cell stage, while at 44 h, almost all 



the uninjected embryos are at the multicell stage. Comparison 
of the earlier timepoint with the later timepoint revealed that 
depletion of a pluripotency factor reduced the extent of zygotic 
genome activation (see Supplementary information and 
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Supplementary Figure S5) . We also found that MO-mediated 
knockdown of Oct4 had the most prominent effect, compared 
with knockdowns of Sall4 and Nanog. At the 20-h timepoint, 
a total of 634 genes were differentially expressed in Oct4 
MO-injected embryos, while <10 genes were differentially 
expressed in Sall4 MO- and Nanog MO-injected embryos (false 
discovery rate (FDR) < 0.1) (Figure 2A) (see Materials and 
methods for details). However, at 44 h post injection, a clear 
effect on global expression levels was observed for all three 
transcription factors (Figure 2B) . Notably, 505 out of 799 genes 



that were downregulated during Sall4 knockdown were also 
downregulated during Oct4 knockdown and 133 out of 245 
genes that were downregulated during Nanog knockdown 
were also downregulated during Oct4 knockdown, indicating 
that the three transcription factors shared a close functional 
relationship during pre-implantation embryonic development 
(P< 10 " chi-square test) . 

We asked which of the differentially expressed genes 
are directly controlled by Oct4, Sall4, or Nanog. Since 
chip experiments could not be readily performed on 
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Figure 2 Microarray analysis of pooled embryos identified the protein-coding regulons of Oct4, Sall4, and Nanog. (A) Venn diagram showing the number of genes that 
were differentially expressed 20 h after MO injection to knock down one of the three pluripotency factors. At this earlier timepoint, only knockdown of Oct4 exerted a 
strong effect on global gene expression levels. (B) Venn diagram showing the number of genes that was differentially expressed 44 h after MO injection. At this later 
timepoint, knockdown of any of the three transcription factors caused hundreds of genes to be significantly downregulated or upregulated (FDR < 0.1). (C) GO analysis 
of genes that were differentially expressed during knockdown of Oct4, Sall4, or Nanog in pre-implantation embryos but showed no evidence of being controlled by the 
pluripotency factors in ESCs. These genes may be novel direct targets of the three transcription factors in the embryos or they may also be indirectly regulated. Upper 
panel: GO analysis of the embryo-specific genes that were downregulated. Lower panel: GO analysis of the embryo-specific genes that were upregulated. All the 
biological processes that are unique to the embryo (P< 0.01) are plotted. The actual numbers of genes for each process (with the expected numbers given in brackets) 
are shown at the end of the bars. Many of the functional categories are related to either transport (in red) or metabolism (in blue). (D) A 15-bp motif within Oct4-bound 
regulatory regions of genes that were also differentially expressed during knockdown of Oct4 in pre-implantation embryos or ESCs. (E) A 7-bp core motif was found to be 
enriched in Nanog-bound regulatory regions of genes that were also differentially expressed during knockdown of Nanog in pre-implantation embryos or ESCs. The core 
motif was extended at the 5' end for genes that were downregulated during Nanog knockdown and it was extended at the 3' end for genes that were upregulated during 
Nanog knockdown. 
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pre-implantation embryos, we compared our microarray 
results with publicly available chIP data sets in ESCs instead 
(Boyer et al 2005; Loh et al 2006; Lim et al 2008; Yang et al 
2008; Chen et al 2008a; Ouyang et aU 2009). 20.4% of the 
downregulated genes and 18.8% of the upregulated genes 
have promoters or enhancers that are bound by the corre- 
sponding transcription factor (Supplementary Figure S6). 
Gene Ontology (GO) analysis using the Panther classification 
system (Thomas et aU 2003} further revealed that the down- 
regulated genes are involved in the estabhshment or main- 
tenance of chromatin architecture, while the upregulated 
genes are involved in stress response and differentiation or 
developmental processes. Many of the direct targets are also 
involved in the cell cycle or are components of key signaling 
pathways (Supplementary Figure S7; Supplementary File S6). 

Since ESCs do not exist in a natural developmental context, 
the pluripotency factors may control additional genes in pre- 
implantation embryos but not in ESCs. To determine the 
embryo-specific genes, we subtracted away from our micro- 
array results genes that were present in the ESC chIP data sets 
as well as genes that responded to Oct4, Sall4, or Nanog 
depletion in ESCs (Supplementary Figure S6}. GO analysis 
revealed that many of these embryo -specific genes were 
associated with metabohsm- and transport-related functions 
(Figure 2C; Supplementary File S7) , which reflected the unique 
embryonic conditions that are absent from ESCs. 

To identify ci5-acting elements that mediate the direct targeting 
of Oct4, Sall4, or Nanog to the genome in pre-implantation 
embryos, we searched the bound loci for known motifs and also 
performed de novo motif discovery using CisGenome (Ji et al 
2008). The Oct4-Sox2 consensus sequence was present in the 
promoter or enhancer regions of 83 downregulated genes (out of 
284 genes) and 35 upregulated genes (out of 153 genes) 
(Figure 2D; Supplementary File S8). In addition, we obtained 
with high confidence a core motif, TGACCTT, that was enriched in 
the Nanog-bound promoter or enhancer regions of 25 down- 
regulated genes (out of 40 genes) and 31 upregulated genes (out of 
46 genes) (Supplementary File S9). Notably, the motif for 
downregulated genes is extended at the 5' end by two basepairs, 
while the motif for upregulated genes is extended at the 3' end 
by one to three basepairs (Figure 2E). A de novo search for 
ci5-regulatory elements using a different program, SCOPE, also 
returned the same core motif (Carlson et al, 2007; Chakravarty 
et al, 2007; Supplementary Figure S8). Taken together, our results 
suggest that different cofactors may mediate Nanog's function as a 
transcriptional activator or repressor. Furthermore, the core motif 
is similar to previously identified motifs for the retinoic acid 
receptor Kara as well as the estrogen-related receptors Esrra and 
Esrrb (Supplementary Figure S9), which suggests that Nanog may 
co-regulate many genes together with Rara and Esrra/Esrrb. 



MicroRNAs form an integral component of the 
Oct4-Sall4-Nanog regulatory network in 
pre-implantation embryos 

MicroRNAs (miRNAs) are an important class of post-tran- 
scriptional regulators that have diverse and critical roles 
during animal development and regeneration (Yang and Wu, 
2007; Stefani and Slack, 2008). It is likely that they also have 
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major roles in pre-implantation developmental processes. We 
therefore quantified miRNA expression from the 1-cell to the 
blastocyst stage in uninjected mouse embryos and detected 
297 distinct miRNAs in at least one of these stages 
(Supplementary File SIO). K-means clustering further revealed 
that many of the miRNAs were temporally regulated over 
development (Figure 3 A) . 

We hypothesize that some of the highly expressed miRNAs 
are integral components of the Oct4-Sall4-Nanog gene 
regulatory network in the pre-implantation embryo. To 
determine which miRNAs are regulated by Oct4, Sall4, or 
Nanog, we knocked down each transcription factor separately 
and measured miRNA expression levels 46 h after microinjec- 
tion when the knockdown phenotypes were readily observa- 
ble. More miRNAs were downregulated than upregulated 
(Figure 3B; Supplementary File Sll), indicating that Oct4, 
Sall4, and Nanog served primarily as positive regulators of 
miRNA expression in pre-implantation embryos. In addition, 
most of the downregulated miRNAs were significant (P< 0.05) 
in at least two different knockdown conditions, suggesting that 
the expression of many miRNAs was under combinatorial 
control by the transcription factors (Figure 3B). In total, 120 
distinct miRNAs were differentially expressed. 40 of them 
(33.3%) have promoters or enhancers that are occupied by 
Oct4, Sall4, or Nanog in ESCs (Lim et al 2008; Marson et al 
2008). The remaining 80 miRNA genes may be regulated 
indirectly by the three transcription factors or they may 
represent novel direct targets of Oct4, Sall4, or Nanog in the 
pre-implantation embryo. 

Next, we sought to understand how the pluripotency factors 
and their protein-coding and non-coding regulons are inter- 
connected. 61.1% of the protein-coding genes (2646 out of 
4329 genes) that were differentially expressed during knock- 
down of Oct4, Sall4, or Nanog were predicted by the miRanda 
algorithm (www.microrna.org) (Betel et al 2008) to be 
targeted by at least one miRNA that was also regulated directly 
or indirectly by the three transcription factors (Supplementary 
File S12) . Similar results were obtained when we restricted our 
target predictions to protein-coding genes and miRNAs that 
were present in the ESC chIP data sets (62.7% ; 622 out of 992 
genes) (Supplementary File SI 3). In comparison, when we 
randomly selected 10 000 distinct sets of protein-coding genes 
and miRNAs, we found that only 52.8 ± 0.7% of the protein- 
coding genes were predicted to be miRNA targets. Hence, our 
analysis indicates that Oct4, Sall4, and Nanog and their 
protein-coding and non-coding regulons are more intercon- 
nected than random (P< 10 " ^, Z-test) . To increase confidence 
in the miRNA-target predictions, we further determined which 
of the miRNAs outputted by miRanda were also predicted by at 
least one other algorithm (see Materials and methods and 
Supplementary Files S14 and SI 5). An example of a subnet- 
work is illustrated in Figure 3C. 



Knockdown of Oct4, Sal 14, or Nanog causes gene 
expression changes even in embryos that appear 
to develop normally 

As our microarray data were generated from pools of embryos 
and the ensemble averages may have masked interesting 
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Figure 3 Oct4, Sall4, and Nanog control the expression of numerous developmentally regulated mlRNAs in the pre-implantation embryo. (A) K-means clustered 
expression profiles for 235 miRNAs that were clearly expressed in uninjected (normal) pre-implantation embryos. Genes are in rows, with temporal progression from left 
to right. Red indicates high expression levels, while green indicates low expression levels. The majority of these miRNAs show a general increase in expression levels 
(cluster I and cluster V) from the 1 -cell to the blastocyst stage. (B) Venn diagram showing the number of miRNAs that were differentially expressed during knockdown of 
Oct4, Sall4, or Nanog (P<0.05). MiRNAs in the miR-290-295 cluster were differentially expressed under all three knockdown conditions. (C) An example of an 
integrated network controlled by Oct4, Sall4, and Nanog in pre-implantation embryos. The protein-coding genes are depicted in various colors to signify different 
functional categories. Dark brown: pluripotency transcription factors; light blue: epigenetic regulators; light brown: signaling molecules; pink: differentiation factors; 
purple: embryo-specific factors. The listed miRNAs are within the top 20 most abundant miRNAs in the embryos. X ^ Y means 'X regulates Y.' For the arrows emanating 
from the center, a solid arrow indicates confirmed direct regulation, while a dotted arrow indicates either indirect regulation or potential direct regulation. 



embryo-to-embryo variability or important expression 
dynamics at the single embryo level, we sought to monitor 
gene expression levels in individual embryos using BioMark 



(Fluidigm). We focused our single embryo analysis on 48 
protein-coding genes, 43 of which were found from our 
microarray experiments to be differentially expressed during 
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Oct4, Sall4, or Nanog knockdown. The 43 genes were chosen 
because they perform essential functions in ESCs, are key 
epigenetic regulators or signaling molecules, have important 
roles during differentiation events, or may have embryo- 
specific functions. The five remaining genes had unchanged 
expression levels and served as negative controls. 

In our first set of analyses, we compared between the 
various experimental conditions, namely (1) no injection, 
(2) injection with a control MO, and (3) injection with a MO 
targeting Oct4, Sall4, or Nanog. For each condition, we 
collected between 80 and 120 embryos at 44 h after MO 
injection. When all the embryos were considered in the 
analysis, we found that our BioMark and microarray results 
agreed reasonably well with each other (Supplementary 
Information, Supplementary Figure SIO, and Supplementary 
File S16}. However, since depletion of a pluripotency factor 
caused embryos to arrest at different developmental stages 
(Supplementary Figure Sll), some of the transcript changes 
observed between experimental conditions could be due to 
expression changes between the various stages. Hence, we 
next compared 4-cell and multicell control embryos with 
knockdown embryos that were also at the 4-cell and multicell 
stages (red boxes in Figure 4A}. Intriguingly, we found that 
many genes were still significantly downregulated or upregu- 
lated (Figure 4B). In particular, at least four genes, namely 
Jarid2, Aqp3, Atg2a, and Sf3b2, were differentially expressed 
in all three knockdown conditions. Hence, although some of 
the knockdown embryos appeared to be morphologically 
normal, numerous molecular changes had already occurred in 
them. 



Expression of DnmtSb distinguishes early- 
arrested embryos from late-arrested embryos 
during Oct4, Sall4, or Nanog knockdown 

Strikingly, we observed that, while uninjected and control- 
injected embryos were able to develop in a stereotypical 
manner, embryos that lacked any of the three pluripotency 
factors arrested over a range of developmental stages instead 
(Supplementary Figure Sll). At 44 h after MO injection, the 
knockdown embryos were distributed mainly over the 2-cell, 
3 -cell, 4-cell, and more advanced multicell stages. Hence, in 
our second set of analyses, we asked if our BioMark 
experiments could identify genes that would differentiate 
early-arrested embryos (2-cell and 3-cell embryos) from late- 
arrested embryos (4-cell and multicell embryos) during 
knockdown of a pluripotency factor (blue boxes in Figure 4A) . 

Interestingly, the expression of a handful of genes correlates 
with the extent to which a knockdown embryo develops. 
Figure 4C shows a heatmap of the top 10 genes that are 
differentially expressed between early-arrested embryos and 
late-arrested embryos. Most of them encode DNA-binding 
proteins, such as the transcription factor Tbx3, which shares 
many common targets with Oct4 and Nanog (Han et aU 2010) . 
We further compared 4-cell knockdown embryos with multi- 
cell knockdown embryos (Table I) and found that 7 of the 10 
genes that were differentially expressed between early- 
arrested embryos and late-arrested embryos were also more 
highly expressed in multicell embryos than in 4-cell embryos. 
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Hence, the expression of these seven genes (highlighted in 
green in Figure 4C and Table I) is low in 2-cell and 3-cell 
knockdown embryos, high in 4-cell knockdown embryos, and 
even higher in multicell knockdown embryos. 

Since promoter DNA methylation contributes to ESC gene 
regulation and methylation-deficient ESCs are restricted in 
their developmental potential (Fouse et aZ, 2008), we decided 
to investigate the DNA methyltransferase Dnmt3b in greater 
depth. Knockdown of Dnmt3b using a translation-blocking MO 
revealed a more severe phenotype than that for any of the three 
pluripotency factors. Specifically, 66.3 ±21.6% of embryos 
injected with the Dnmt3b MO arrested by the 2-cell stage, even 
though the percentages for Oct4, Sall4, and Nanog were only 
between 5.5 and 19.0% (P<0.001, Student's f-test) (Figure 5A). 
Rescue experiments with in vitro transcribed Dnmt3b mRNA 
confirmed the specificity of the Dnmt3b MO (Figure 5B-D). 
These results suggest that Dnmt3b is essential for pre- 
implantation development and that it may have a role in 
determining the extent to which an embryo depleted of Oct4, 
Sall4, or Nanog can develop. In addition, although Dnmt3b 
embryos had been reported to develop normally before 
E9.5 (Okano et al 1999b; Ueda et al 2006), it is likely that 
maternal Dnmt3b, which we detected in 1-cell mouse zygotes 
(Supplementary Figure SI), may be driving the earliest 
developmental stages in the knockout mouse model. 



A coherent feed-forward loop controls the 
expression of DnmtSb 

In ESCs, chip experiments showed that the pluripotency 
factors can control the transcription of Dnmt3b directly. 
Interestingly, the factors can also regulate Dnmt3b expression 
indirectly via the ESC-specific miRNAs in the miR-290-295 
cluster. Specifically, Oct4 and Nanog directly control the 
expression of these miRNAs, which shows a general upward 
trend from the 1-cell to the blastocyst stage (Supplementary 
Figure S12). The miRNAs then target the 3'UTR of the 
retinoblastoma-like 2 (Rbl2) mRNA, whose protein product 
functions as a transcriptional repressor of Dnmt3b (Benetti 
et aU 2008; Sinkkonen et aU 2008). In the pre-implantation 
embryo, we found from our microarray analysis that although 
Rbl2 did not pass our cutoff criterion for significant genes, its 
expression is still upregulated 1.6- to 2.9-fold during knock- 
down of Oct4, Sall4, or Nanog. Our microarray analysis also 
revealed that Dnmt3b expression was downregulated by 3- to 
19-fold when Oct4, Sall4, or Nanog was depleted (FDR<0.1; 
Supplementary Files S3 and S5). Hence, a coherent feed- 
forward loop (FFL) consisting of the pluripotency factors, 
the miR-290-295 cluster, Rbl2, and Dnmt3b regulates 
DNA methylation in the pre-specified embryo and in ESCs 
(Figure 8B). 

To obtain a quantitative understanding of the feed-forward 
architecture, we constructed a model of non-linear ordinary 
differential equations (ODEs) based on Michaehs-Menten 
kinetics with cooperativity (Figure 6A) and performed the 
simulation using a Runge-Kutta algorithm. We chose the 
parameters so that the model would reproduce known 
experimental observations. In uninjected (wild-type) 
embryos, the pluripotency factors are present at high levels 
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and activate the expression of the miR-290-295 cluster. The 
miRNAs then downregulate the production of Rbl2. When the 
level of Rbl2 is low, Oct4, Sall4, and Nanog can then fully 
activate the expression of Dnmt3b (Figure 6B, left panel). 
However, the converse is true during knockdown of the 
pluripotency factors (Figure 6B, right panel). Furthermore, in 
Dicer-null ESCs, both the transcript and protein levels of 
Dnmt3b were significantly lower than that in wild-type ESCs 
(Benetti et al, 2008; Sinkkonen et al, 2008) . In agreement with 
these observations, when we set the expression of the miRNAs 
to zero in our model, we found that the peak Dnmt3b 
expression level dropped by 1.7-fold for our chosen set of 
parameters (Supplementary File SI 7). 

We next investigated the advantages that the FFL confers on 
the pluripotency factors' regulation of DNA methylation. 
Superficially, the pathway through the miRNAs and Rbl2 
appears to be unnecessary, as the pluripotency factors can 
simply activate the expression of Dnmt3b. To identify 
differences between this straightforward design and the FFL, 
we also simulated the two-component model consisting of the 
pluripotency factors and Dnmt3b only (Supplementary File 
SI 8). When we varied the concentration of the pluripotency 
factors from low (0) to high (1) and examined the output 
(normalized Dnmt3b expression level) of both models, we 
found that the simple two-component model produced a more 
graded output, while the FFL gave a more switch-like response 
(Figure 6C). Hence, the miRNA-Rbl2 pathway served to 
sharpen the system output so that it has a more switch-like 
behavior. 

Since the performance of a non-linear dynamical system is 
dependent on the precise values of the parameters, we 
wondered how our results may change with different 
parameter choices. We first found that the feed-forward 
network motif is robust with respect to the synthesis rate of 
the Rbl2 transcriptional repressor (see Supplementary 
information and Supplementary Figure S13). Next, we studied 
the behavior of the dynamical system when we varied the 
synthesis rate of Dnmt3b. The expression of Dnmt3b is under 
combinatorial regulation by the pluripotency factors and Rbl2. 
In our mathematical description of the FFL, we have modeled 
the influences of each transcription factor as both mutually 
independent (second and third terms on the right-hand side of 
the third ODE) and cooperative (fourth term of the third ODE) . 
However, the relative contributions of each regulator to the 
transcription of Dnmt3b are unknown. Thus, we varied the 
ratio of the two parameters, (3a and and observed how the 
concentration of Dnmt3b responds to changes in the level of 
the pluripotency factors. We found that when Pa/Pb is 



approximately bigger than 5, the hypersensitive response is 
lost and the FFL essentially behaves like the simple two 
component system. However, when the ratio becomes smaller, 
the switch-like behavior becomes more pronounced 
(Figure 6D; Supplementary Figure SI 4). In contrast, the 
steepest slope of the response curve is unaffected by the 
parameter Pc (see Supplementary information and 
Supplementary Figure SI 5). Hence, our analysis showed that 
in order to achieve a hypersensitive response, the regulation of 
Dnmt3b transcription cannot be dominated by the pluripo- 
tency factors, so that Rbl2 can exert its effect. 

A prediction from our model is that at high levels of 
pluripotency factors, for example in pre-specified embryos 
and in ESCs, the concentration of Dnmt3b becomes more 
susceptible to intrinsic gene expression noise when the FFL is 
broken. To test the prediction, we determined the cell-to-cell 
variability of Dnmt3b transcript level for wild-type ESCs and 
for Dicer"/" ESCs (Calabrese et al, 2007). Specifically, 
we performed quantitative real-time PCR (RT-qPCR) experi- 
ments on four different sets of cells: (1) Dicer"/" ESCs, (2) 
Dicer"/" ESCs transfected with a control miRNA mimic that 
has minimal sequence identity with miRNAs in humans and 
mice (Dharmacon), (3) Dicer " / " ESCs transfected with a miR- 
295 mimic (Dharmacon), and (4) wild-type ESCs. To verify 
that the transfection of miR-295 was successful, we measured 
the levels of its target genes Rbl2, p21, and Casp2 (Benetti etal, 
2008; Sinkkonen et al, 2008; Wang et al, 2008; Zheng et al, 
2011) and confirmed that they were all downregulated 
(Supplementary Figure SI 6). We also transfected each of the 
four cell populations with a plasmid encoding green fluor- 
escent protein (GFP), so that we could sort for single cells by 
flow cytometry, and then determined the transcript levels of 
Dnmt3b, Oct4, Sall4, and Nanog in individual ESCs using 
Biomark (Fluidigm) . Strikingly, Dnmt3b expression exhibited 
the widest spread of Ct values in Dicer-null ESCs. The 
interquartile range for Dicer-null ESCs and for cells transfected 
with the control miRNA was 8.94 and 3.33 respectively, 
while the range for Dicer-null ESCs transfected with miR-295 
and for wild-type ESCs was 2.87 and 2.79 respectively 
(Figure 6E, upper panel) . The standard deviation of Dnmt3b 
expression level was also larger in Dicer-null ESCs and control- 
transfected cells than in miR-295-transfected cells and wild- 
type ESCs (Supplementary Figure S17A). Furthermore, these 
results are not due to the expression of the pluripotency factors 
exhibiting different degrees of variability between the four sets 
of cells (Figure 6E, lower panel, and Supplementary Figure 
S17B). When corrected for the standard deviation of the 
pluripotency factors, the expression of Dnmt3b was still more 



Figure 4 Gene expression analysis of individual embryos reveals hidden molecular changes and transcript differences between early- and late-arrested embryos. 
(A) A schematic diagram showing the developmental progression of uninjected and injected embryos. Embryos were collected at 44 h post-injection for expression 
analysis. Two main analyses were performed. First, we compared the expression of control embryos that were at the 4-cell and multicell stages with the expression of 
knockdown embryos that were also at the 4-cell and multicell stages (red boxes). Second, we compared the expression of knockdown embryos that were at the 2-cell and 
3-cell stages (early-arrested) with the expression of knockdown embryos that were at the 4-cell and multicell stages (late-arrested) (blue boxes). (B) The first comparison. 
Although all the embryos in this analysis were morphologically similar (4-cell and multicell stages), many genes were still differentially expressed in the knockdown 
embryos. The individual Ct values for the top 10 ranking genes of each knockdown condition are plotted as heatmaps. P-values are given in brackets. Top panel: Oct4 
knockdown; middle panel: Sall4 knockdown; bottom panel: Nanog knockdown. Most of the genes shown are downregulated upon knockdown of a pluripotency factor. 
The genes that are upregulated are marked by an asterisk (*). In addition, genes that are highlighted in pink are common among all the three knockdown conditions and 
they participate in critical functions such as epigenetic gene silencing (Peng et al, 2009; Li et al, 2010; Pasini et al, 2010), autophagy (Tsukamoto et al, 2008), and 
alternative splicing (Tang etal, 2010). (C) The second comparison. The individual Ct values for the top 10 ranking genes are plotted as a heatmap, with the P-values 
given in brackets. The expression levels of the genes highlighted in green are also significantly different between the 4-cell stage and the multicell stage. 
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variable in Dicer-null cells and control-transfected cells 
than in miR-295-transfected cells and wild-type ESCs 
(Supplementary Figure S17C}. Hence, our results indicate 

A 



that the FFL serves to buffer Dnmt3b expression against 
intrinsic variability in the concentrations of the pluripotency 
factors. 
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Single-cell analysis reveals an increase in gene 
expression noise during knockdown of Oct4 or 
Sall4 

Our simulations suggest an intriguing explanation for why 
during knockdown of a pluripotency factor, some embryos 
might express a higher level of Dnmt3b and progress further in 



Table I The top 10 genes that were differentially expressed between 4-cell 
knockdown embryos and multicell knockdown embryos 



Gene 


P-value 


Fold change 


Fgfr2 


1.80 X 10"^^ 


2.44 


Hmgbl 


2.13 X 10"^^ 


2.20 


PouSfl 


8.16 X 10"^^ 


2.46 


Dnmt3b 


1.47 X 10"^^ 


1.91 


Tbx3 


1.52 X 10"^^ 


2.45 


Sf3b2 


6.92 X 10"^ 


1.48 


Atg7 


4.53 X 10"^ 


2.02 


Racgapl 


6.17 X 10"^ 


1.49 


Klf5 


8.43 X 10"^ 


1.64 


Dnmt3a 


2.31 X 10"^ 


1.46 



The genes highlighted in green were also expressed at significantly lower levels 
in 2-cell and 3 -cell embryos when compared with later stage embryos. 



development, while other embryos might express a lower level 
of Dnmt3b and arrest earher in development. In uninjected 
embryos, the pluripotency factors are present at high levels 
and fluctuations in their concentrations due to intrinsic noise 
do not affect the level of Dnmt3b significantly (right side of 
blue graph in Figure 6C} . However, when a pluripotency factor 
has been knocked down to a low level, the gradient of the 
response curve increases dramatically and the system exists in 
the switch regime instead (left side of blue graph in Figure 6C} . 
In this case, even small fluctuations in the concentrations of 
the pluripotency factors can affect the level of Dnmt3b 
appreciably. 

To test the hypothesis that the expression of Dnmt3b and 
possibly many other genes become noisier during knockdown 
of a pluripotency factor, we performed single-cell real-time 
PGR experiments using BioMark (Fluidigm) . In these experi- 
ments, we focused on Oct4 and Sall4, since knockdown of 
Nanog gave a milder phenotype that is consistent with the 
expression of Nanog being significantly lower than that of 
Oct4 or Sall4 in the earUest developmental stages. For each 
experimental condition (no injection, injection with control 
MO, injection with Oct4 MO, and injection with Sall4 MO), we 
collected embryos that were morphologically at the 4-cell stage 
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Figure 5 The DNA methyltransferase DnnntSb is essential for pre-implantation development. (A) Percentages of embryos that arrested by the 2-cell stage for four 
different knockdown conditions. Unexpectedly, we observed that the majority of embryos injected with a MO targeting DnmtSb arrested around the maternal-zygotic 
transition (P< 0.001, Student's f-test, when comparing DnmtSb with any of the pluripotency factors). (B) Representative microscopy images of pre-implantation mouse 
embryos under various experimental conditions. (C) Quantification of observed phenotypes after 4 days of /n vWro culture. Compared with embryos that were injected 
with 0.6 mM DnmtSb MO only, co-injection of 0.6 mM DnmtSb MO with 1-10ng/|il DnmtSb mRNA reduced the fraction of embryos that arrested by the 2-cell stage 
(7 independent experiments) (P<0.05, Student's f-test). (D) Percentages of embryos that developed until at least the morula stage after 4 days of culture. Compared 
with embryos that were injected with DnmtSb MO only, co-injection of DnmtSb MO with /n wfro transcribed DnmtSb mRNA increased the percentage of embryos that 
reached the morula stage (P<0.05, Student's f-test). 
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Time (arbitrary units) Time (arbitrary units) 




Normalized concentration of pluripotency factors Ctrl miR miR-295 

Figure 6 A coherent feed-forward motif sharpens the DnnntSb transcriptional output into a bistable switch. (A) A model of non-linear ODEs describing the behavior of the 
FFL. In the model, both the rate of transcription and the miRNAs' silencing of Rbl2 are modeled on Michaelis-Menten kinetics with cooperativity. In the equations, B represents 
basal promoter activity, y stands for the rate of degradation, /C stands for the Michaelis constant, and H is the Hill coefficient. A list of parameter values that reproduces known 
experimental observations is provided in Supplementary File S1 7. (B) Simulations of the system showing how the concentrations of the miRNAs (blue), DnmtSb (red), and Rbl2 
(green) change over time. Left panel: high levels of the pluripotency factors; right panel: low levels of the pluripotency factors. (C) A plot showing how the concentration of 
DnmtSb varies as a function of the concentration of the pluripotency factors. Both the original FFL (blue curve) and a simplified two-component system comprising only the 
pluripotency factors and DnmtSb (red curve) were simulated. Only the FFL exhibits a threshold or switch-like behavior. (D) A plot showing how the DnmtSb transcriptional 
output depends on the relative regulatory strengths of different transcription factors. A larger pa/Pb ''a^io indicates that DnmtSb transcription is more reliant on the pluripotency 
factors, while a smaller pa/pb ratio indicates that the expression of DnmtSb is more dependent on the Rbl2 transcriptional repressor. (E) RT-qPCR results for single Dicer 
knockout (Dcr KO) ESCs, Dcr KO ESCs transfected with a control miRNA mimic, Dcr KO ESCs transfected with a miR-295 mimic, and wild-type ESCs. Upper panel: A boxplot 
for the Ct values of DnmtSb. Lower panel: A boxplot for the Ct values of the pluripotency factors. The bottom and top of every orange box represent the 25th and 75th 
percentiles and the band inside each box is the median. The dotted lines extend to the furthest datapoint that is within 1 .5 times the interquartile range from the box, while the 
black circles indicate the outliers. The Ct values of DnmtSb for the Dicer-null ESCs are significantly larger than those for the miR-295 transfected cells or the wild-type ESCs 
(P< 0.005, Wilcoxon rank sum test) but are comparable to the Ct values for the control-transfected cells (P> 0.01 , Wilcoxon rank sum test). In contrast, the Ct values of the 
pluripotency factors are similar between the Dicer-null ESCs and the other three cell populations (P>0.01, Wilcoxon rank sum test). 
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and manually separated them into individual cells, so that we 
could profile gene expression in each blastomere. All the cells 
were harvested at 39 h post-injection, as uninjected embryos 
were at the 4-cell stage at this timepoint. 

We monitored the expression of 84 functionally diverse 
genes, including Dnmt3b (Supplementary File S19). 68 of 
these genes were selected because they were differentially 
expressed during knockdown of one of the pluripotency 
factors, while the remaining genes were chosen because they 
were related family members or performed similar functions. 
Hierarchical clustering of our single-cell data revealed that all 
the blastomeres from uninjected embryos grouped together in 
a single tight cluster. However, the cells from Oct4 or Sall4 MO- 
injected embryos did not appear to form any coherent cluster. 
We further calculated the standard deviations of Ct values for 
every gene within each experimental condition and then 
plotted a histogram for these standard deviations (Figure 7 A) . 
Notably, the distributions for Oct4- and Sall4-depleted 
blastomeres were shifted to the right, indicating that gene 
expression was more variable in the knockdown cells 
compared with the control cells. The average standard 
deviation for uninjected cells was 1.00, while the average 
standard deviations for Oct4-depleted and Sall4-depleted cells 
were 1.74 and 3.47 respectively and were significantly higher 
than that for uninjected cells (P<10~^ and P<10"^^ 
respectively, Wilcoxon rank sum test) . In contrast, cells from 
control MO-injected embryos showed an average standard 
deviation of only 1.21 (P>0.1, Wilcoxon rank sum test). 
Hence, knockdown of Oct4 or Sall4 resulted in noisier gene 
expression, not just for Dnmt3b but also for the majority of the 
genes that we profiled. 

To confirm the results, we performed principal component 
analysis (PCA) on our single-cell data. The first principal 
component accounted for 70% of the gene expression 
variability between all the cells, while the second principal 
component accounted for 7% of the variability. A projection of 
the cells' expression vectors onto these two principal 
components showed that the uninjected blastomeres and the 
control MO-injected cells clustered together, while the knock- 
down cells were clearly more spread out (Figure 7B}. In 
addition, we observed that among all the knockdown 
blastomeres, cells from different embryos were generally 
further apart from one another, while cells from the same 
embryo were closer to one another, indicating that the gene 
expression from embryo to embryo is noisier than the gene 
expression between cells of the same embryo. Nevertheless, 
increases in both inter-embryo variability and intra-embryo 
variability contributed to noisier gene expression during 
depletion of a pluripotency factor. 

The increase in gene expression noise is not a result of 
variable knockdown efficiency among the different embryos, 
as determined by quantification of protein levels in our 
immuno staining images of control and knockdown embryos 
(Figure lA; Table II; also see Supplementary information and 
Supplementary Figure S18}. We further co-injected each MO 
with the corresponding mRNA and showed that the increase in 
gene expression noise could be rescued (see Supplementary 
information and Supplementary Figure SI 9). Collectively, our 
data demonstrate that consistent depletion of a pluripotency 
factor in preimplantation mouse embryos leads to an increased 
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variability in the expression of numerous genes, including 
Dnmt3b. 



Discussion 

How multiple cellular events are coordinated in the correct 
spatiotemporal order to ensure proper development of the 
1-cell fertilized zygote into a pluripotent blastocyst that can 
give rise to all cell types in the body is a largely unsolved 
problem in mammahan developmental biology. While gene 
regulatory networks have been extensively studied in ESCs in 
the past decade, the limiting amount of material available for 
experiments has hampered the apphcation of chIP to map 
transcription factor binding sites in pre-implantation embryos. 
Furthermore, the technical challenge of microinjecting the 
small and fragile embryos poses a barrier to extensive gene 
knockdown studies. 

Here, we leveraged on recent advances in gene expression 
profiling methods to deconstruct the pluripotency gene 
regulatory network that controls pre-implantation develop- 
ment in the mouse embryo. We knocked down Oct4, Sall4, and 
Nanog by microinjecting MOs into the embryos and profiled 
the gene expression of pooled embryos, single embryos, and 
individual blastomeres to identify the regulon of the pluripo- 
tency factors. Similar to their roles in ESCs, Oct4, Sall4, and 
Nanog control the expression of numerous genes involved in 
the cell cycle, establishment or maintenance of chromatin 
architecture, signaling, and differentiation processes in the 
pre-implantation embryo. We also discovered that many of the 
genes that were differentially expressed in the embryo during 
knockdown of Oct4, Sall4, or Nanog were not regulated 
directly or indirectly by the three pluripotency factors in ESCs. 
GO analysis of these embryo-specific genes revealed that they 
largely performed metabohsm- or transport-related functions. 
Interestingly, a recent study found that the expression of 
metabohsm genes undergoes major changes during the 
derivation of ESCs from the ICM, which may be a result of 
the cells accommodating a dramatic change in their environ- 
ment (Tang et aU 2010} . Hence, since the embryonic environ- 
ment is distinct from the in vitro culture conditions of ESCs, 
Oct4, Sall4, and Nanog may have an extended role in the pre- 
implantation embryo by regulating the unique metabolic 
requirements and cellular conditions of the embryo. 

There are other pluripotency factors that may perform 
important roles in pre-implantation embryos, in particular 
Sox2, which is known to co-regulate many genes with Oct4 
and Nanog in ESCs. However, we found that maternal Sox2 
protein is difficult to deplete. In embryos injected with a MO 
targeting the translation start site of Sox2, we could still detect 
by immuno cytochemistry a substantial amount of Sox2 
protein 24 h after injection. Unsurprisingly, these embryos 
did not show an obvious phenotype until the morula 
stage. Future advances in chIP technologies will allow us to 
study the Sox2 regulatory network in the earhest develop- 
mental stages. 

Strikingly, we observed that embryos injected with Oct4, 
Sall4, or Nanog MOs arrested over a range of developmental 
stages. Uninjected wild-type mouse embryos showed a 
relatively synchronous rate of pre-implantation development 
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Figure 7 Single-cell analysis reveals that gene expression becomes noisier during knockdown of a pluripotency factor. (A) A plot showing how the standard deviations 
of gene expression are distributed for each experimental condition. Gray bars: uninjected embryos, orange bars: embryos injected with a control MO, red bars: embryos 
injected with a MO targeting Oct4, and green bars: embryos injected with a MO targeting Sall4. Most genes showed an increased variability in expression when Oct4 or 
Sall4 was depleted (P< 10" and P< 10" respectively, Wilcoxon rank sum test). We also noted that the gene expression variability of our uninjected and control- 
injected cells is similar to the variability observed in a previously published single-cell profiling data set for wild-type blastomeres (Tang etal, 201 1) (see Supplementary 
information). (B) Principal component projections of individual blastomeres showing how similar or different the cells are to one another. Blastomeres that were derived 
from uninjected 4-cell embryos (solid blue circles) or from embryos injected with a control MO (solid red circles) clustered close together. Blastomeres that were derived 
from embryos injected with a MO targeting Oct4 (unfilled black circles) or Sall4 (unfilled green squares) were more spread out. Each knockdown cell is also labeled as 
'E_B_', where 'E' stands for embryo and 'B' stands for blastomere. From the plot, we observed that embryo-to-embryo variability is noticeably larger than the variability 
between cells from the same embryo. For example, for the Sall4 knockdown condition, E2B1 , E2B2, E2B3, and E2B4 are positioned relatively close to one another but 
they are all far away from any blastomere from E3 (also see Supplementary information). 



Table II Variability in protein level (standard deviation/mean) based on 
quantification of pixel intensities 



Antibody 


Uninjected 


Control injected 


Knockdown 


Background 


a-0ct4 


0.172 


0.123 


0.166 


0.167 


a-Sall4 


0.170 


0.284 


0.292 


0.082 



with over 90% developing to the blastocyst stage by Day 4 of 
in vitro culture. However, the majority of embryos injected 
with Oct4, Sall4, or Nanog MOs arrested at 2-cell, 3-cell, 4-cell, 
and multicell stages. This loss in developmental control is not 
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due merely to the act of injecting and possibly injuring the 
embryos because, even though embryos injected with a 
control MO developed shghtly slower than uninjected 
embryos, over 80% of them still progressed to at least the 
morula stage by Day 4. 

Why did some knockdown embryos arrest early, while 
others arrest at a later stage? From our single-cell analysis, we 
discovered that depletion of a pluripotency factor caused an 
increase in gene expression noise. The increase is contributed 
predominantly by larger variations in gene expression 
between different embryos and also partly by a larger spread 
of expression levels between individual cells of the same 
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embryo. As a result, many genes are expressed at high levels in 
some embryos but at low levels in other embryos. A subset of 
these genes might be driver genes that perform key functions 
and determine the developmental potential of an embryo. 
Hence, embryos that express higher levels of the driver genes 
will progress further in development, while embryos that 
express lower levels of the driver genes will arrest earher in 
development (Figure 8A} . 

We found that the DNA methyltransferase Dnmt3b may 
function as one of the driver genes. The expression of Dnmt3b 
correlates strongly with the extent to which an embryo 
depleted of Oct4, Sall4, or Nanog can develop. To determine 
whether it is a driver gene, we separately knocked down 
Dnmt3b and discovered that over 60% of the knockdown 
embryos arrested by the 2-cell stage, which is around the time 
when the embryonic genome becomes activated in the mouse. 
In addition, a recent study showed that, although the pre- 
implantation embryonic genome is in a general state of 
hypomethylation, the oocyte actually contributes a substantial 
number of CpG island promoters that remain hypermethylated 
in the pre-specified embryo (Smith et al, 2012). Hence, it is 
possible that Dnmt3b is required to maintain the methylation 
status of these regions during and after the maternal-to- 
embryonic transition. 

A FFL regulates the expression of Dnmt3b in pre-implanta- 
tion embryos and in ESCs, whereby pluripotency factors can 
control the expression of Dnmt3b directly or they can regulate 
its expression indirectly through the miR-290-295 cluster of 
miRNAs and the Rbl2 transcriptional repressor (Figure 8B}. 
Modeling and simulation of the FFL provided insights into the 
underlying cause of the increase in gene expression noise 
during knockdown of a pluripotency factor. Compared with a 
straightforward two-component design where the pluripo- 
tency factor simply activates Dnmt3b in the absence of other 
pathways, the feed-forward motif allows the expression of 
Dnmt3b to exhibit a more switch-like behavior, so that the 
system approaches bistability. Importantly, we note that the 
switch between an 'off state and an 'on' state occurs when the 
expression level of the pluripotency factors is low. For our set 
of parameter values (Supplementary File SI 7), simulation 
showed that when the expression of the pluripotency factors is 
high, a 20% fluctuation in their levels causes a 3% change in 
the expression of Dnmt3b. However, when the expression of 
the pluripotency factors is low, a 20% fluctuation in their 
levels causes a 60% change in the expression of Dnmt3b 
(Figure 6C}. Hence, in an uninjected embryo where the levels 
of Oct4, Sall4, and Nanog are high, the transcriptional output 
of Dnmt3b is robust to intrinsic noise in the expression of its 
activators. However, in knockdown embryos, the transcrip- 
tional output of Dnmt3b becomes much more sensitive to 
small changes in the concentration of its activators. 

While the FFL may serve as an important paradigm for how 
pre-implantation mouse embryos ensure their relatively 
synchronous rate of development, we speculate that there 
are many other undiscovered network features in the embryo 
that either buffer its development against intrinsic and 
extrinsic noise or to provide safeguard measures. Like the 
miR-290-295 cluster, we found that the most highly expressed 
miRNA in the pre-implantation embryo, miR-709, also 
depends on Oct4 and Sall4 for its expression. Interestingly, a 
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recent study showed that miR-709 targets and degrades BORIS 
mRNA, which has been imphcated in the genome-wide 
erasure of DNA methylation in the germline (Tamminga 
et aU 2008). Hence, the high expression of miR-709 may 
provide insurance against the accidental erasure of essential 
methylation marks in the embryo (Figure 8B) . 

Two closely related de novo DNA methyltransferases, 
Dnmt3a and Dnmt3b, are encoded in the mouse genome. 
Here, we have focused on Dnmt3b for a few reasons. 
First, although the expression of both Dnmt3a and Dnmt3b 
was assayed in our single-embryo experiments, only Dnmt3b 
allowed us to distinguish early-arrested embryos from 
late-arrested embryos. The expression of Dnmt3b was 2.1- 
fold higher in embryos that progressed to the 4-cell and 
multicell stages than in embryos that arrested at the 2-cell and 
3-cell stages (Figure 4C), while the expression of Dnmt3a was 
only 1.1-fold higher (P> 0.001). Second, Dnmt3a~^" mice 
developed to term and appeared grossly normal at birth. In 
contrast, the Dnmt3b-null allele resulted in embryonic 
lethality, indicating a more prominent role for Dnmt3b during 
embryogenesis (Okano et aU 1999a). Third, there is only one 
known promoter for the Dnmt3b gene, but there are two 
separate promoters for the Dnmt3a gene (Chen et aU 2002), 
which can comphcate analysis. For example, the level of the 
full-length Dnmt3a protein is similar between Dicer ESCs 
and wild-type ESCs, while the level of the smaller Dnmt3a 
protein, which originates from a downstream intronic promo- 
ter, is appreciably lower in Dicer-null ESCs (Sinkkonen et aZ, 
2008). Clearly, the regulatory activities at the two promoters 
can be different and further work is needed to tease apart the 
regulation occurring at the two distinct Dnmt3a promoters. In 
the future, chIP assays to map the genome-wide locations of 
the DNA methyltransferases and gene knockdown experi- 
ments followed by bisulfite sequencing will enable us to 
understand how Dnmt3a and Dnmt3b contribute to shaping 
the pre-implantation methylome. 

In conclusion, we have deconstructed an Oct4-Sall4-Nanog 
gene regulatory network in pre-implantation mouse embryos 
that coordinates various critical cellular events to ensure 
rehable developmental progression despite intrinsic gene 
expression noise. Not only have we delineated functional sets 
of genes dependent on the pluripotency factors for their 
expression, we have also demonstrated how a particular 
network design, the FFL, may buffer embryos against 
fluctuations in the concentration of the hub proteins, so that 
the embryos can develop in a relatively synchronous and 
stereotypical manner. Our work has laid the foundation for 
future studies on other gene networks controlling pre- 
implantation mammahan embryogenesis and also has impor- 
tant imphcations for assisted reproduction. 

Materials and methods 
Embryo collection and culture 

All procedures involving animals were performed under our Institu- 
tional Animal Care and Use Committee (lACUC) protocol #8315 titled 
'Novel roles of pluripotency regulators in the early mouse embryo', 
which was approved by the Administrative Panel on Laboratory 
Animal Care (APLAC) at Stanford University. Three- to five-week-old 
wild-type Fl [C57BL6 x DBA/2) females (Charles River) were 
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Figure 8 Schematics of an exemplary regulatory circuit governing developmental progression in the pre-implantation embryo and the consequences of its 
perturbations. (A) A graphical model outlining the sequence of events that occurs during knockdown of Oct4, Sall4, or Nanog. DnmtSb is a representative driver gene, 
whose expression level may determine the extent to which a pre-implantation embryo can develop. (B) An Oct4-Sall4-Nanog regulatory network controlling DNA 
methylation in the pre-implantation embryo. The FFL, which is one of the most prevalent motifs or basic building blocks of transcriptional networks (Shen-Orr etal, 2002; 
Milo et al, 2002; Mangan and Alon, 2003; Shalgi et al, 2007; Tsang et al, 2007; Gerstein et al, 2012), is highlighted by a yellow box. 



superovulated with intraperitoneal injections of 5 lU pregnant mare's 
serum gonadotropin (PMSG; Sigma) followed by 5 lU human chorionic 
gonadotropin (hCG; Sigma) 48 h later before being mated overnight 
with wild-type Fl males. One-cell zygotes were released from oviducts 
17 h after hCG injection and cumulus cells were removed by 
hyaluronidase (Sigma) treatment and pipetting. Embryos were 
cultured in drops of Human Tubal Fluid with 10 % serum supplement 
(SAGE In- Vitro Fertilization Inc.) under mineral oil (Sigma) in 5 % CO2 
at 37°C. Each 20 \x\ drop contained 8-10 embryos. 



Single-cell isolation 

The zona pellucida was first removed by acid tyrode's solution and 
the embryos were then incubated in TrypLE solution (Invitrogen) for 
5-10 min at 37°C. The dissociation of individual cells from whole 
embryos was monitored under the microscope. After rinsing twice in 
drops of PBS containing 3 mg/ml polyvinylpyrrolidone (PBS-PVP), the 
dissociated blastomeres were carefully transferred into sample 
collection tubules with the aid of a finely pulled glass pipette. 



Microinjection of antisense MOs 

25-nt antisense MOs were purchased from Gene Tools, LLC. The 
sequences for the MOs are as follows: 

Oct4-M01: 5'-AGT CTG AAG CCA GGT GTC GAG CCAT-3' 
Oct4-M02: 5'-CTC CGATTT GCATAT CTG GGC AGG G-3' 
Sall4-M01: 5'-TCC GCC ACC AAT TCC TGG ACT TG-3' 
Sall4-M02: 5'-AAG AGC CCT GGT GAG TTG CCC CCT C-3' 
Nanog-MOl: 5'-GGA CCA GGA AGA CCC ACA CTG ATG T-3' 
Nanog-M02: 5'-GCA ACA ACC AAA AAA CTG ACT GTC T-3' 
Dnmt3b-M0: 5'-GAT GTC TGC TGT CTG CCT TCATTG T-3' 
Control-MO: 5'-TCC AGG TCC CCC GCATCC CGG ATC C-3' 
The control MO does not target any known sequence in the mouse 
genome or transcriptome. For each control or knockdown condition, 
5-10 pi of 0.6 mM MO was injected into the cytoplasm of each zygote 
on an inverted microscope (Olympus 1X70) equipped with a hydraulic 
micromanipulation system (Narishige IM300 microinjector). In every 



experiment, at least 8-10 embryos were used for each of the 
conditions, namely uninjected, injection of control MO, or injection 
of MO targeting one of the pluripotency factors. 



Immunocytochemical staining 

Embryos were first fixed in 4% paraformaldehyde for 30 min and then 
permeabilized in 0.4% Triton X-100 for 30 min. Next, they were 
blocked in 10 % normal goat serum for 1 h before incubation with 
antibodies. The embryos were incubated in 1:100 diluted primary 
antibody overnight at 4°C, in 1:5000 diluted secondary antibody for 
1 h, and in 3 |iM DAPI for 10 min, before being mounted on slides using 
Fluoromount-G (SouthernBiotech). All antibodies were diluted in 
blocking solution. The primary antibodies used were mouse mono- 
clonal anti-0ct4 (Santa Cruz Biotechnologies, sc-5279), mouse 
monoclonal anti-Sall4 (Santa Cruz Biotechnologies, sc- 101147), and 
rabbit polyclonal anti-Nanog (Cosmo Bio, RCABOOOIP). Secondary 
antibodies were purchased from Invitrogen/Molecular Probes (Alexa 
Fluor 594 donkey anti-mouse IgG (A-21203) and Alexa Fluor 594 
donkey anti-rabbit IgG (A-21207)). Embryos were washed in between 
the various steps with PBS-PVP and were imaged on the Leica SP2 
AOBS confocal laser scanning microscope. 



Synthesis of mRNAs for rescue experiments 

The overall strategy was to clone the coding region from plasmids 
using primers that contained the T7 promoter. The PGR products then 
served as templates for in vitro transcription, which was performed 
using the mMessage mMachine T7 Ultra Kit (Ambion). 

For Oct4, the plasmid MMM1013-98478834 (Open Biosystems) was 
first modified so that the final mRNA cannot be targeted by the 
corresponding MOl. The third basepairs of the DNA sequence coding 
for the second to eighth amino acids were changed using the 
QuikChange Lightning Kit (Agilent) and the following mutagenic 
primers: 5'-CCC ACC TTC CCC ATG GCA GGT CAT GTC GCATCT GAT 
TTC GCC TTC TCA CCC CC-3' and 5'-GGG GGT GAG AAG GCG AAA 
TCA GAT GCG AGA TGA CCT GCC ATG GGG AAG GTG GG-3'. The 
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template for the transcription reaction was then amplified from the 
modified plasmid using the following primers: 5'-CCA AG T AAT ACG 
ACT CAC TAT AGG GA C AAG AGG TAATCC ACC CCC GGC TCG-3' and 
5'-CAA AAT GAT GAG TGA GAG ACA GGC GAG GCT-3' (the T7 
promoter is underlined) . 

For Sall4, the plasmid OMM5896-99848745 (Open Biosystems) was 
used for cloning. Even though the gene insert does not include the MO 
target region, it lacks a STOP codon. The STOP codon and a FLAG tag 
were inserted at the C'-terminus via two rounds of site-directed 
mutagenesis using the QuikChange Lightning Kit (Agilent). The 
primers for the first round of mutagenesis were 5'-TAA GAT TGC 
TGT GAG CGA CGA CAA GTA ATC AGG CCT CAT GGG CC-3' and 
5'-GGC CCA TGA GGC CTG ATT ACT TGT CGT CGC TGA CAG CAA TGT 
TA-3', while the primers for the second round of mutagenesis were 
5'-AAATAA GAT TGC TGT CAG CGA CTA CAA GGA CGA CGA CGA CAA 
GTA ATC AGG CC-3' and 5'-GGC CTG ATT ACT TGT CGT CGT CGT CCT 
TGT AGT CGC TGA CAG CAA TGT TAT TT-3'. The template for the 
transcription reaction was then amplified from the modified plasmid 
using the following primers: 5'-CCA AG T AAT ACG ACT CAC TAT AGG 
GCG CGC ACC ATG TCG AGG CGC AAG CAG GC-3' and 5'-AAA GCT 
GGG CCC ATG AGG CCT GAT TAC TT-3'. 

For Nanog, the plasmid EMM1002-99866141 (Open Biosystems) was 
modified so that the final mRNA cannot be targeted by 
the corresponding MOl. The third basepairs of the DNA sequence 
coding for the second to seventh amino acids were changed using the 
QuikChange Lightning Kit (Agilent) and the following mutagenic 
primers: 5'-CCATCA CAC TGA CAT GAG CGT CGG ACT ACC AGG ACC 
CCA CAG TTT GCG TAG T-3' and 5'-ACT AGG CAA ACT GTG GGG TCG 
TGG TAG TCG GAC GCT CAT GTC AGT GTG ATG G-3'. The template for 
the transcription reaction was then amplified from the modified 
plasmid using the following primers: 5'-CCA AG T AAT ACG ACT CAC 
TATAGGGTT GCG TAA AAC CTT TTC AGA AAT CCC TTC C-3' and 
5'-CCG ACT GCT CTT CGG AAG GTC AGG AGT TC-3'. 

For Dnmt3b, the plasmid MMM1013-99827219 (Open Biosystems) 
was also modified so that the final mRNA cannot be targeted 
by the corresponding MO. The third basepairs of the DNA sequence 
coding for the second to seventh amino acids were changed using the 
QuikChange Lightning Kit (Agilent) and the following mutagenic 
primers: 5'-CCA GCG TCA CGA CAG GAA ACA ATG AAA GGC GAT AGT 
AGG CAC CTG AAT GAA GAA GAG GGT G-3' and 5'-CAC CCT CTT CTT 
CAT TCA GGT GCG TAC TAT CGC CTT TCA TTG TTT CCT GTC GTG 
AGG CTG G-3'. The template for the transcription reaction was then 
amplified from the modified plasmid using the following primers: 
5'-CCA AG T AAT ACG ACT CAC TAT AGG G CG GCG CAA GTA AAC GTA 
GCG CAG CGAT-3' and 5'-GGG GAA CGA TGT GCT TGT GGT CCT AGC 
TTC-3'. 



Cell culture and transfections 

The Dicer + (wild-type) and Dicer"/" ESCs (Calabrese a/, 2007) 
were grown in knockout DMEM (Life Technologies) supplemented 
with 15% FBS (Thermo Scientific Hyclone), Glutamax (Life Technol- 
ogies), non-essential amino acids (Life Technologies), penicillin/ 
streptomycin (Life Technologies), p-mercaptomethanol (Life Technol- 
ogies) , and LIF (Millipore) . Cells were transfected with Lipofectamine 
2000 (Life Technologies) according to manufacturer's instructions. 
The media was replaced 4-6 h after transfection. 



Microarray analysis of pooled embryos 

For each experiment, RNA was extracted from 20 embryos using the 
PicoPure RNA Isolation Kit (Arcturus). The RNA was reverse 
transcribed and amplified using the WT-Ovation One-Direct Amplifi- 
cation System (NuGEN). The cDNA was then fragmented and labeled 
using the FL-Ovation cDNA Biotin Module (NuGEN) and hybridized 
to the Mouse Genome 430 2.0 Array (Affymetrix) . All manufacturers' 
instructions were followed. The expression data were submitted to 
the NCBl Gene Expression Omnibus (GEO) (http://www.ncbi. 
nlm.nih.gov/geo) under accession number GSE42135. 

The raw data of 42 GeneChips were normalized and processed using 
dChip (Li and Wong, 2001). The estimations of expression level, fold 



change between two samples, and the FDR were calculated based on 
the methodology developed by Johnson and Wong (see http:// 
www.stanford.edu/group/wonglab/doc/NewGeneScoreBiostatDraft. 
pdf). In comparison with the traditional moderated t-statistic, this 
methodology takes into account both the fold change and the 
expression magnitude and calculates a statistic based on fold change, 
logged fold change, and unpaired t-statistic. We defined differential 
expression by a threshold of 0.1 FDR estimated from 5000 random 
permutations across the samples. At least one of the following two 
conditions has to be satisfied for a gene to be considered as being 
regulated by a pluripotency factor: (1) significant in uninjected sample 
versus knockdown sample but not significant in uninjected sample 
versus control injected sample or (2) significant in control injected 
sample versus knockdown sample. We converted accession IDs of the 
results to gene names using tools available on the Mouse Genome 
Informatics (MGI) website. 



Analysis of chIP sequencing (chlP-seq) data sets 

The genomic regions bound by Oct4 and Nanog were obtained from a 
data set published by Chen et al (2008b) . The genes regulated by Oct4 
and Nanog were obtained in two ways. First, we searched for the 
closest gene to the bound chromosomal position. Second, we obtained 
a list of regulated genes using the approach developed by Ouyang et al 
(2009) with a transcription factor-association strength (TEAS) scores 
threshold of 3. The list of genes directly regulated by Sall4 was pooled 
from two published data sets (Lim et al, 2008; Yang et al, 2008). The 
accession IDs from Yang et al were converted to gene names using the 
tools on the MGI website. We used CisGenome (Ji et al, 2008) to 
calculate motif enrichment and to perform de novo motif discovery 
based on mouse genome mm8. The bound loci from Chen et al were 
given as single nucleotides, so we extended them to 400 bp regions 
centered at those loci (200 bp upstream and 200 bp downstream) for 
our motif analysis. Control regions were randomly selected and the 
number of control regions was five times that of the targets. 



GO analysis 

The Gene Expression Data Analysis tool within the PANTHER resource 
(http://www.pantherdb.org/tools/genexAnalysis.jsp) (Thomas et al, 
2003) was used to identify functional categories in our microarray 
results. The default set of Mus Musculus genes was used as the 
reference. Since some common categories are 'noisy' and will always 
appear regardless of the input list, we also randomly generated subsets 
of genes from all the genes that were tiled on the Affymetrix microarray 
and performed GO analysis on these random gene lists. The categories 
or biological processes that showed up repeatedly during the GO 
analysis of these random lists were then discarded. 



Expression profiling of mlRNAs 

For each experiment and each condition, 20 embryos were harvested 
directly into DNasel-containing Lysis Solution from the TaqMan 
MicroRNA Cells-to-Ct Kit (Applied Biosystems) and the lysates were 
stored at - 80°C. We used the Megaplex™ Primer Pools (Rodents; 
Applied Biosystems) for reverse transcription and pre-amplification. 
Reverse transcription was performed at 42°C for 30min using the 
MultiScribe RTand pre-amplification was performed for either 12 or 16 
cycles (95°C for 15 s and 60°C for 4min) using the TaqMan PreAmp 
Master Mix (Applied Biosystems). Real-time PGR was performed on 
the 7900HT System (Applied Biosystems) using the TaqMan Rodent 
MicroRNA Array Set (also called TaqMan Low Density Arrays; Applied 
Biosystems). 

Our pipeline for converting Ct values to expression levels is as 
follows. We set a threshold Ct of 32, above which we considered the 
miRNA to be absent. The data were normalized to the average Ct 
values for the U6 RNA and the expression level was taken to be 2 " 
For each biological sample, we repeated the experiment twice, once 
with 12 cycles of pre-amplification and once with 16 cycles of pre- 
amplification. If the expression level from 12 cycles of pre-amplifica- 
tion was 0, then we kept the expression level from 16 cycles of 
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pre-amplification unless there was evidence of non-specific amplifica- 
tion (typically an unusually small Ct outlier). Otherwise, 
we would use the expression level from 12 cycles of pre-amplification, 
as there would be less amplification bias with fewer cycles. To 
compare knockdown embryos with uninjected embryos or embryos 
injected with a control MO, we performed a regularized f-test (Baldi 
and Long, 2001) to determine the differentially expressed miRNAs 
(P<0.05). 



MiRNA-target predictions 

We first downloaded a set of target site predictions from www.mi- 
crorna.org (August 2010 release), which was based on the miRanda 
algorithm (Betel et al, 2008), to identify the miRNAs that may 
potentially regulate the protein-coding genes that were differentially 
expressed during knockdown of a pluripotency factor. Subsequently, 
to increase confidence in the predictions, we used the miRWalk 
database (Dweep et al, 2011) to identify the miRNA-target pairs that 
were predicted not only by miRanda but also by at least one other 
program, including DIANA-microT (Maragkakis et al, 2009a, b), 
miRDB (Wang and El Naqa, 2008; Wang, 2008) , miRWalk (Dweep et al, 
2011), PITA (Kertesz et al, 2007), RNA22 (Miranda et al, 2006), and 
TargetScan (Lewis et al, 2005; Friedman et al, 2009). Default settings 
were used for the minimum seed length and the P-value. 



IVIicrofluidic quantitative real-time PCR 

For reverse transcription and pre-amplification, TaqMan gene expres- 
sion assays (20 x , Applied Biosystems) were pooled, so that each 
assay was at a final concentration of 0.2 x . Each individual embryo 
was harvested directly into a RT-PreAmp Master Mix (5.0 |il CellsDirect 
2 X Reaction Mix (Invitrogen) , 2.5 |il 0.2 x Assay Mix, 0.5 |il Super- 
Script 111 RT/Platinum Taq Mix (Invitrogen), 1 .4 |il TE Buffer, and 0.1 |il 
SUPERase-ln (Ambion)). The harvested samples were first reverse 
transcribed at 50°C for 20min and then pre-amplified for 20 cycles 
(95°C for 15 s and 60°C for 4min). The pre-amplified products were 
diluted four-fold with TE buffer before real-time PCR on the BioMark 
System (Fluidigm). We used 48.48 Dynamic Arrays (Fluidigm), 
TaqMan Universal PCR Master Mix (Applied Biosystems), and TaqMan 
gene expression assays (Applied Biosystems) in these real-time PCR 
experiments. The arrays were primed, loaded, and run according to 
manufacturer's instructions. Ct values were calculated from the 
BioMark System's software using default settings. 



Single-embryo and single-cell data processing 

To test for normality, we applied the Shapiro-Wilk test. If the two data 
sets under consideration (either original Ct values or relative 
expression levels) passed the test, then we would use the Student's 
t-test to compare them. If not (which means either one or both data sets 
did not come from a normal distribution), then we would apply the 
Wilcoxon Rank Sum test. In addition, we applied a two-sample 
Kolmogorov-Smirnov test for all comparisons. Expression levels were 
calculated using a AACt method. 

Since individual cells are derived from whole embryos, fluctuations 
in single-cell gene expression values can be due to embryo-to-embryo 
variability (inter-embryo variability) or due to differences between 
cells of the same embryo (intra-embryo variability) . To determine the 
contribution of inter-embryo variability to gene expression fluctua- 
tions observed between individual cells, we first computed the total 
expression level of an embryo from the Ct values of the four cells that 
constituted that embryo and then calculated the standard deviation 
between whole embryos. To determine the contribution of intra- 
embryo variability to single-cell gene expression fluctuations, we first 
computed for each gene the standard deviation of the set of Ct values 
from cells of the same embryo and then calculated the mean of these 
standard deviations for each experimental condition (uninjected 
embryos, embryos injected with a control MO, or embryos injected 
with a MO targeting a pluripotency factor) . 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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